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Abstract 

High-resolution neutron backscattering techniques are exploited to study the elastic and quasi- 
elastic response of the high-density amorphous (HDA), the low-density amorphous (LDA) and 
the crystalline ice I c upon temperature changes. Within the temperature ranges of their structural 
stability (HDA at T < 80 K, LDA at T < 135 K, ice I c at T < 200 K) the Debye-Waller factors and 
mean-square displacements characterise all states as harmonic solids. During the transformations 
HDA^LDA (T rj 100 K), LDA^I C (T « 150 K) and the supposed glass transition with T g an 135 K 
no relaxation processes can be detected on a time scale t < 4 ns. It can be concluded from coherent 
scattering measurements (D2O) that LDA starts to recrystallise into ice I c at T ~ 135 K, i.e. at 
the supposed T g . In the framework of the Debye model of harmonic solids HDA reveals the highest 
Debye temperature among the studied ice phases, which is in full agreement with the lowest Debye 
level in the generalised density of states derived from time-of-flight neutron scattering experiments. 
The elastic results at low T indicate the presence of an excess of modes in HDA, which do not obey 
the Bose statistics. 

PACS numbers: 63.50. +x, 64.70.Kb, 78.70.Ck 
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I. INTRODUCTION 



Water, although comprising one of the most simple molecules, reveals extremely com- 
plicated properties in its condensed states. In the solid state for example, there exist not 
only twelve well known crystalline phases but also a number of disordered structures 
Two of the disordered structures are known as high-density amorphous (HDA) and low- 
density amorphous (LDA) ice, characterised by a molecular density of p ~ 39 molec./nm 3 
and p ~ 31 molec./nm 3 , respectively . The origin of a number of different disordered 
phases within a single substance, often referred to as amorphous polymorphism, has not 
been understood, yet. It has been conjectured on the basis of computer simulation re- 
sults that HDA and LDA might represent the supercooled glassy structures of two different 
liquid phases of water Obviously, this scenario implies the existence of a 

glass-transition designated by a glass-transition temperature T g . 

From the experimental point of view, there is a number of results backing the findings from 
computer simulations. Structures equivalent to HDA and LDA can be obtained from the 
liquid by extremely fast quenching of /zm-droplets or from the gas phase by vapour deposition 
Q y H Q. As far as LDA-type structures are concerned, they all transform 
exothermally to a crystalline cubic phase (ice I c ) at T rs 150 — 160 K on an endothermal 



plateau setting in at T « 135 K, which has been o 



(DSC) and other calorimatric experiments 
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orimetry 



26|. The 



very onset of the endothermal contribution has been interpreted as water's T g in literature, 
whereby in LDA the lowest T g « 124 K of all low-density amorphous structures is discussed 

Q. 

A comparable T„ has been deduced from DSC studies of different aqueous solutions, 



whose specific T g s extrapolate upon dilution toward the single value of 135 K 
Finally, from observation of D 2 0-H 2 interlayer mixing upon heating of LDA-type vapour 



deposited water a higher molecular mobility at T « 150 K has been concluded [2i 

However, these experimental results do not establish unequivocally a glass-transition 
or a glass-transition temperature in water. For each of the cited experiments there can 
be found at least one experiment reported in literature whose results either contrast the 
glass transition scenario or give an alternative explanation for the observed phenomena. For 
example, endothermal transitions have been observed at the same thermodynamic conditions 
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in transformations of high-density crystalline ice phases to ice I c . They are understood as 
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entropy driven order — > disorder transitions in the proton sublattice of the phases |25 . 
a feature equally detected in ice In this context, it is noteworthy that recently well 



established DSC results on apparent LDA samples 
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23, 



24j had to be reinterpreted 



as such a reversible order — > disorder transition in the proton sublattice of ice XII 



. The 



very low T g « 124 K reported in remains to be reexamined 2^|. Consequently, one may 
ask what is the unique dynamic feature by which, on the one hand, the endothermic plateau 
in the amorphous systems is interpreted as a glass-transition but, on the other hand, it is 
identifiied in the crystalline structures as a disorder phenomenon of the proton sublattice 
only. Whether the recently observed endothermic effect in pressure dependent experiments 
on emulsified amorphous samples is a unique feature of the disordered system requires equally 
a closer examination 26^ . 

In contrast to the D 2 H 2 interlayer mixing effect 
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30 1 ion diffusion experiments have 

n 

not shown any higher mobility of guest molecules in the amorphous ice matrix 34] neither 
lave isotopic exchange measurements shown any translational motion of water molecules 
. Instead, they have given the picture of a defect promoted proton mobility responsible 
for the endothermic transition, a feature directly detectable by Nuclear Magnetic Resonance 

Recently, models of a so called shadow glass-transition have been proposed that explain 
the thermal anomalies observed at T « 135 K with a real but experimentally unaccessible 



glass-transition at T E « 165 K ySj. Whether this shadow glass-transition could be indeed 

n n 

applicable to water is the issue of current debates [39, 40]. 

Another idiosyncratic feature which has not been satisfactorily explained, yet, is the in- 
elastic response of the amorphous ice structures. Dynamic properties of bulk HDA, LDA 
and LDA-type structures differ strongly from established features of supercooled glassy sys- 



tems 0,|42j. The dynamic response of LDA does not exhibit any phonon dampinc 



any remarkable Boson-peak 



4J,|45J,|4 



£0], and any two-level systems (TLS) |48l. l4fl 



44]. 



Consequently, its spectral density and thermal conductivity are reminiscent of a harmonic 
crystalline state (51 1. In the case of HDA, this departure from the characteristics of glassy 
states is apparently less pronounced, since TLS modes have been unequivocally observed 



in optical absorption experiments 



neutron experiments 



HQ 



48, 



49 



50]. 



Similar features are observed by inelastic 



, where an excess of modes in the density of vibrational 
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states is reported. However, the detailed phonon dispersion and also the spectral density at 
T > 40 K of HDA are identified as crystal-like in other experiments pointing at an intriguing 



high degree of short-range order 43J, |55 ] . 



In the light of the above outlined experimental facts the idea of HDA and LDA being 
strongly disordered crystalline systems not thermodynamically connected to water's liquid 



phase |56J], may at first g 
been discussed for water 



ance seem rather tempting. Such a nano-crystalline scenario has 



561 ] but also for other systems showing amorphous polymorphism 



. A discontinuity between the liquid state, vapour deposited and hyper-quenched amor- 
phous water on the one hand and LDA on the other hand has been indeed postulated based 
on theoretical concepts and results from spectroscopic experiments and computer simulations 



58, 
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Giving the subject some further thinking it becomes, however, obvious that 
a nano-crystalline scenario does not provide a more stringent explanation for the absence 
of glassy features in the dynamics. The inherently high disorder within the nano-crystals 
should leave a clear trace in the low-frequency dynamics not too different from that of a 
glassy state. 

To advance the question further it is obviously indispensable to characterize the funda- 
mental, microscopic relaxation processes responsible for the endothermic plateau in DSC 
experiments. This characterization must be carried out both in time and space. While 
order-disorder phenomena and the glass relaxation may take place on similar time scales 
the spatial correlations are different. In particular, the spatial patterns should allow to 
discern the hypothetical relaxation channels of the glass-transition from processes observed 
during the transitions of high-density crystalline structures to ice I c . 

The present work represents our first attempt of scrutinizing the glass-transition hypoth- 
esis of water with a supposed T g « 135 K on a microscopic scale by neutron backscattering 
and time-of-flight techniques, i.e. by spectroscopic techniques sampling time and spatial 
correlations. Both techniques have been extensively used for studying glass transition phe- 
nomena and proved to be substantial for their understanding. Both, time-of-flight and 
backscattering spectroscopy offer the opportunity of studying the microscopic dynamics of a 
sample directly as an energy resolved response. In addition, neutron backscattering can be 
also used to exploit the elastic scattering within a narrow energy range ~ 1 /ieV from which 
important dynamic properties can be deduced |d^ . 0|. This enables us, in principle, to 
characterise relaxation processes on a molecular scale during the apparent glass-transition. 
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All experiments performed reveal that the properties of HDA and LDA follow a tem- 
perature dependence in complete accordance with the harmonic theory of the solid state 
63) 3- Beyond the low-temperature limit, where zero-point oscillations are predominant, 
the molecular mean-square displacements of HDA and LDA are proportional to the temper- 
ature increase. This behaviour is expected if the degrees of freedom, and thus the spectral 
density of the sample in its different states is determined by harmonic modes only The 
very good conformity with the mean-square displacements deduced from the phonon den- 
sities of states of HDA and LDA reinforces this scenario. In the temperature regime of the 
transitions HDA — > LDA and LDA — > ice I c comprising the apparent glass-transition with 
T g ~ 135 K no higher mobility of water molecules, e.g., in terms of translational diffusion, is 
observed on a time scale shorter than some nano-seconds. From experiments on coherently 
scattering samples the temperature T « 135 K is identified as the onset of a recrystallisation 
of LDA into ice I c jf^. In addition, an excess of modes in HDA, that however does not 
follow the Bose-statistics valid for harmonic vibrations, is clearly indicated in the values of 
the zero-point oscillation. 

To bring these findings clearly forward the present paper is structured in the following 
way. The next two chapters render some useful details on the applied neutron scattering 
techniques, the samples and experimental procedure, and introduce some observables 
(Debye- Waller factor, mean-square displacement, density of states, velocity of sound), to 
which we refer throughout the text. In section IV, we present and discuss the experimental 
data in view of the studied transitions (HDA — > LDA, LDA — > I c , and the glass-transition), 
the Debye- Waller factor and mean-square displacement, and the excess of modes in HDA. 
All results obtained are summarised in section V. 



II. EXPERIMENTAL 



To meet the requirement of high energy resolution the experiments have been performed 
on the neutron backscattering spectrometers IN13 and IN16 and the time-of-flight spec- 
trometer IN6 at the Institut Laue Langevin in Grenoble, France jsg. The principles of the 
backscattering spectrometers are based on the neutron beam monochromatization and the 
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energy analysis of the scattered beam by single crystal Bragg reflection 



n-A = 2dsin(6) (f) 



in backscattering geometry O rs 90° 



62 



63j. Using the CaF 2 (422) reflection on IN13 



and the Si(lll) on IN16 the incident neutron energies (E) and energy resolutions (AE) of 
E = 16.5 meV, AE = 8/ieV and E = 2 meV, AE = 1/xeV are obtained, respectively. Beyond 
the differences in incident energy and energy resolution these spectrometers are optimised 
for sampling complementary Q-ranges. IN16 with Q < 2.1 A" 1 is best utilised for e.g. long 
range diffusion processes with a correlation length of r > 3 A. IN13 with Q < 5.5A -1 is 
optimised for short range diffusion with r > 1 A. Thus, with both instruments we are able 
to cover a spatial range stretching from correlations within the proton sublattice of water 
up to its intermolecular distances. 

In general, backscattering spectrometers can be utilised in two different modes. Firstly, in 
the elastic scan mode, the incident neutron energy is kept constant to permanently meet the 
Bragg condition of the analyser given by eq. The elastic scan mode samples changes of the 
elastic intensity within AE due to dynamic and relaxation processes which may be induced 
by changing experimental parameters like temperature in the present case. Secondly, in the 
energy scan mode, the incident neutron energy is varied in a well defined, systematic way 
within a narrow energy range 8E, e.g., 8E ~ ±10yueV on IN16. In the energy scan mode, 
dynamic and relaxation processes can be, in principle, characterised quantitatively. Both, 
energy and elastic scan mode, were applied on the spectrometer IN16 whereas IN13 was 
used in the elastic scan mode exclusively. 

The time-of-flight spectrometer IN6 was used with incident energies of 3.1 meV (AE = 
80 fjeV) and 4.8 meV (AE = 150 /ieV). The sampled Q-ranges of 0.3-2.1 A" 1 (3.1 meV) 
and 0.3-2.6 A -1 (4.8 meV) matched the Q regime of the IN16 measurements. However, the 
energies probed by the instrument IN6 correspond to dynamic processes on a time scale of 
t < 50 pico-seconds. 

In detail, we studied two fully protonated samples on IN16, one sample, S3, in elastic 
scan and a second, S4, in energy scan mode, and two samples on IN13, sample S2 fully 
deuterated and sample SI partially protonated with 40% vol. H 2 0. Three fully deuterated 
samples were measured on IN6. The deuteration was employed for a better control of the 
sample state via the pronounced coherent scattering contribution of D 2 to the signal and 
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for an enhanced signal contribution from the Oxygene atoms. Figure ^ shows the elastic 
intensity I(Q, ou — 0) of the samples SI and S2 in the HDA and the LDA state in comparison 
to the structure factor S(Q) = f S(Q,u)du of fully deuterated HDA and LDA measured 
with high resolution on the diffractometer D20 at the Institut Laue Langevin j^2|. The 
coherent signal of the partially protonated sample is, as expected, strongly suppressed and 
the scattering characteristics are reminiscent of the incoherent contribution. 

All samples were prepared by compression of crystalline ice Ih at T « 77 K (in liquid N2) 
in a piston-cylinder apparatus up to p ~ 18 kbar 6sj|. The formed HDA was recovered at 
ambient pressure, filled into Aluminum holders used as standard in neutron scattering and 
placed into precooled (T = 75 K) standard cryostats. After remaining N 2 had been carefully 
removed at T = 78 K the samples were cooled down to T = 2 K. Table E gives details on 
the thermal treatment and thermal cycling applied to the samples during the experiments. 
Data were accumulated on IN13 for 60 min per point with the samples gradually heated 
corresponding with a rate of 3.5 K/60 min. IN16 elastic scan data were accumulated for 2 min 
with a heating rate of 1 K/2 min. IN16 inelastic scan data were accumulated for 2 x 90 min 
with a heating step of 2 K in the range of the supposed glass transition, whereby in situ 
data updates were performed every 10 min giving, within the data statistics, identical results 
as obtained after 90 min acquisition time. As for the LDA state the thermal treatments 
resulted in phases which had been annealed for 2 h at 125 K plus 1 h at 128.5 K on IN13 and 
for 1 h at 130 K and a preceeding slow heating on IN16 before the regime of the apparent 
glass transition was reached. Thus, annealing times and heating rates corresponded with 
treatments performed in DSC experiments on water's T g reported in literature. 

Samples measured at IN6 were subject to a comparable thermal cycling procedure 
whereby the LDA structure was annealed at 127 K and ice I c was annealed at 160 K. Mea- 
surements following the glass-transition have been performed at the temperatures T = 127 K 
for 4 x 60 min., at T = 132 K for 4 x 60 min., at T = 137 K for 5 x 30 min., at T = 140 K for 
6 x 30 min., at T = 143 K for 4 x 30 min., at T = 147 K for 30 min. and at T = 155 K for 
30 min. before annealing the formed ice I c at 160 K. Long time measurements of 4-6 hours 
were carried out with HDA at 2 K, 20 K, 40 K, 60 K and 80 K, with LDA at 2 K, 20 K, 
40 K, 60 K, 80 K and 127 K, and ice I c at 127 K and 160 K. 

Standard data corrections were performed, accounting for empty container and back- 
ground scattering, self-shielding, absorption effects of the samples and for efficiencies of 
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TABLE I: Thermal cycling applied to the samples SI and S2 (IN13) and S3 and S4 (IN16). The 
phases given below indicate the states in which the samples existed when heat treatment was 
started at 2 K for the measurements. Please note that SI was measured twice within the stability 
range of HDA. 

phase sample SI sample S2 sample S3 sample S4 

HDA 2^65^2^125 K 2^125 K 2^130 K 2^160 K 

LDA 2^170 K 2^170 K 2^260 K 

ice I c 2^170 K 



the instrument detectors. However, due to the scattering power of the samples (optimised 
for multiple scattering behaviour) and the negligible absorption of the scatterers only mi- 
nor perturbations of the signal were observed. Nevertheless, the elastic signal measured in 
the low-Q range on the spectrometer IN16 encounters a lower sensitivity to temperature 
changes in comparison to the results taken on IN13. This feature is well reported in litera- 
ture and associated with contribution from multiple scattering processes 69]. Depending on 
the information which is to be extracted from the experimental data, different normalisation 
standards have to be used. These standards are explicitely described in the text below. 

For a clear presentation of the extensive amount of data some data sets have been re- 
grouped to give statistical errors of the order of the data symbols. Error bars are plotted 
with the data in the figures. The grouping is explicitly indicated in the figure captions for 
data whose interpretation exceeds a pure qualitativ description. All fitting routines were 
performed with ungrouped data sets. The simplex minimization algorithm was applied 
having taken into account the x-grid and y-error of the fitted data. 

Please note, that we focus primarily on the backscattering measurements due to the 
superior resolution of the instruments, when discussing potential relaxation processes 
during the phase transitions HDA to LDA and LDA to ice I c . 



III. OBSERVABLES IN NEUTRON SCATTERING 



For the convenience of the reader let us summarise some basic principles of neutron 
scattering as they are needed for the comprehension of the data analysis and the discussion. 
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FIG. 1: Bottom: Elastic intensity I(Q, uj = 0) of the amorphous samples measured on the spectrom- 
eter IN13 at T ~ 2 K. Solid and dotted lines represent the HDA and LDA structures studied with 
the deuterated sample S2, respectively. ■ and □ represent the HDA and LDA structures studied 
with the partially protonated sample SI, respectively. Top: For comparison static structure factor 
of the deuterated HDA and LDA samples measured with high resolution on the diffractometer 
D20. 

A more detailed introduction into this field can be found in textbooks on neutron scattering 



|3 63, Q 



63l l64j . At the beginning we restrict our consideration to a harmonic system whose 
scattering characteristic is assumed to be isotropic. These terms are well met by amorphous 
solids. For such a system the intensity measured in an elastic scan mode neutron scattering 
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experiment can be written as 



S(Q, u = 0) = exp (W(Q, T)) ■ S(Q) (2) 
W(Q,T) = -Q 2 (u 2 (T)) (3) 

with W(Q,T) the Debye- Waller factor, (u 2 (T)) the mean-square displacement of the scat- 
terers, and S(Q) the static structure factor, which may be approximated for incoherent 
scattering by a scaling factor (2ttN) ~ 1 representing the number of scatterers in the sample 
Harmonicity implies that (u 2 (T)) can be related to the density of states denoted as 
G(u) as 

WT)) = -- /^cothf-^V- . (4) 

M representing the mass of the scatterers jzi|- Taking the low and high temperature limits 
into consideration (u 2 (T)) can be asymptotically approximated by the moments of energy 
(uj^ 1 ) und (uj~ 2 ) as 

<^°» = £7t^ • (5) 

As a consequence, the mean-square displacement is expected to cross over from a constant 
value at low temperatures to a oc T behaviour upon heating leading to the relation 

W(Q,T) oc -Q 2 ■ T . (7) 

Please note, that (u 2 (T — > 0)) measured in neutron scattering experiments is not the absolute 
zero-point oscillation of the scatterers. Its value depends on the incident neutron energy. 
This energy sets a limit to the energy range on which the spectral density of the scatterers 
is sampled at T — ► K. However, (u 2 (T — > 0)) denotes a well defined quantity which makes 
possible to draw comparative conclusions upon the properties of different sample states 
under identical experimental conditions. 

Going back to eq. 121 we can quantitatively determine two physical effects which can induce 
changes in the observable elastic intensity S(Q,u = 0). Firstly, the structure of the sample 
detectable as S(Q) can change due to phase transformations. Secondly, the spectral density 
G(uS) can alter and influence both the low and high temperature properties of S(Q,u) via 
eqs. El and respectively. From the experimental point of view neutron scattering offers 
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the opportunity of discerning these two points by utilising, on the one hand, incoherent 
scatterers (protonated samples SI, S3 and S4) to study changes in G(u) and, on the other 
hand, coherent scatterers (deuterated samples S2) to measure also changes in S(Q). 
Finally, utilising the Debye approximation for G(u) [3] 

with v the average velocity of sound 



,12 

TTTT ' ~3 (9) 



v 3 vf vf 



the observables G{uS) and (u 2 (T)) given in eqs. El and |H1 can be compared qualitatively 
and quantitatively with results from other experimental techniques accessing the velocity of 
sound. Within the Debye model a characteristic quantity is the Debye temperature Td 

with 3sN/V the density of vibrational modes of the system, which can be also extracted by 
numerical techniques from eq. |U as 

o fc2 T 2 r T D /T 

W T )> = 2Mk^'TE' J x - coth ( x / 2 )- dx ■ ( n ) 



IV. RESULTS AND DISCUSSION 

IV. a. Phase transformations and glass— transition 

To give an overview on the temperature dependence of the elastic signal and its changes 
during the thermal treatment, Fig. |21 reports on the intensity I(u = 0) measured in elastic 
scans with SI on IN13 and S3 on IN16. I(u = 0) is the intensity integrated over the 
accessible Q range of the instruments. The data sets are normalised to the intensity detected 
at T = 2 K in the crystalline cubic phase I c (SI, IN13) and LDA (S3, IN16). The phase 
transformations HDA — > LDA and LDA — > I c can be identified in both samples as rather 
sharp intensity gains in the elastic response at T « 110 K and T w 150 K. They are 
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indicated by vertical arrows. Additionally, a small step in intensity at T « 240 K indicates 
the formation of the hexagonal crystal from ice I c as reported in reference jl^l . 73^ . 

At sufficiently high T (> 40 K), the intensity of all phases decreases linearly upon heating 
in the logarithmic presentation of Fig. |21 This behaviour is stressed by solid lines. Neither 
in the T range of HDA— >LDA nor in the T regime of the LDA— T c and, thus, the assumed 
glass-transition, a pronounced drop off in I{uj = 0) can be detected. Such a drop off would 
be a mandatory fingerprint for an enhancement in molecular mobility of the sample (eq. |2J) 
on time scales less than a few nano-seconds. 

A missing enhancement in molecular mobility is also confirmed by energy scan measure- 
ments with the incoherently scattering sample S4 on IN16. Spectra measured at different 
temperatures are reported in Fig. |H1 To compensate for the here nonrelevant temperature 
dependence of the elastic intensity the data are normalised to the maximum intensity at 
the respective temperatures. The identity of all spectra demonstrates that within the in- 
strumental resolution of 1 /ieV no relaxation processes are present in the sample. From 
the comparison of the different data to the results measured at 2 K it is evident that the 
spectra do not change markedly their shape in the energy window of ±10 /ieV, regardless 
of the state of the sample (The contribution of phonons in this energy window is too weak 
to leave a trace in the temperature dependence.) In accordance with the IN16 data, no 
change of the inelastic response apart from the normal harmonic phonon contribution can 
be observed in the IN6 experiments at temperatures prior to the recrystallisation of the 
amorphous samples. We may conclude that according to eqs. El and El there is no evidence 
for a development of additional decay channels for the elastic intensity at the assumed glass 
transition with T g « 135 K. There is no evidence of relaxation behaviour in the course of 
the transformations HDA — ► LDA and LDA — ► I c , either. Consequently, the presence of 
relaxation processes like translational diffusion over intermolecular distances on a time scale 
shorter than some nano-seconds can be excluded. 

Before we will discuss the Debye- Waller factor in detail and, thus, the potential relaxation 
processes on shorter length scales, it is important to establish the regime of stability of the 
amorphous structures. We focus here primarily on the supposed glass transition and the 
LDA — > I c transformation. Fig. 0] reports on the elastic intensity I(Q,u = 0) measured with 
S2 on IN13 at different temperatures which correspond to the regime of the glass transition. 
The temperatures are indicated by vertical arrows in Figure |U c. Here, I(Q,u = 0) is 
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FIG. 2: Temperature dependence of the elastic intensity I(oj = 0) measured with the partially 
protonated sample SI on IN13 and the fully protonatcd sample S3 on IN16. The signal is integrated 
over the accessible Q-range of the instruments. Data on HDA are represented as (IN13) and 
■ (IN13 and IN16), on LDA as □ (IN13 and IN16), and on cubic ice I c as ♦. Details on the 
thermal treatment of the samples are given in the text. The intensity is normalised to the I(u = 0) 
at T ss 2 K of I c on IN13 and of LDA on IN16. Upward arrows indicate the transformation 
HDA — > LDA at T ~ 110 K, downward arrows indicate the transformation LDA — > I c at T ~ 150 K 
and from cubic to hexagonal ice at T ss 240 K measured on IN16 exclusively. Solid lines stress the 
linear behaviour of the intensity expected for a harmonic solid at sufficiently high temperature. 
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FIG. 3: Energy scans I(Q,lo) of the fully protonated sample S4 measured on IN16. All spectra are 
normalised to the maximum intensity determined at the respective temperatures that are indicated 
in the figures. Top figure reports spectra measured in the course of the HDA — > LDA transfor- 
mation, bottom figure reports spectra measured in the temperature range of the apparent glass 
transition. Also plotted are data obtained at T s=s 2 K representing the instrumental resolution. 
Please note, that there is no difference between the I(Q,u>) at different temperatures. 
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FIG. 4: Left: Temperature evolution of the elastic intensity I(Q, to = 0) during the transformation 
of LDA into ice I c measured on the spectrometer IN13 between T ~ 80 K and T ~ 170 K. The 
nominal temperatures are indicated by arrows in figure c in the right panel. The two arrows 
indicate the Q-values at which the temperature dependence I(T,u = 0) is inspected in detail in 
figures a and b in the right panel. Right: Elastic intensity I(T,oj = 0) measured at Q = 1.7 A -1 
in figure a, at Q = 2.4 A" 1 in fi gure b, and integrated over the accessible Q-range in figure c. 
Dashed vertical line stresses the supposed glass transition and the grey shaded area indicates the 
transition regime of LDA to ice I c 

normalised to I(Q,u = 0) of LDA determined at 2 K. It gives, therefore, the relative 
changes before and after the conversion of LDA into I c . The two arrows indicate the Q 
values at which the temperature dependent intensity I(T,u = 0) is plotted in Figure |U a 
and b. As it is expected from a harmonic system I(Q,u = 0) exhibits up to 130 K only 
a loss of intensity due to the Debye- Waller factor (eq. |HJ). Whereas at 170 K I(Q,u = 0) 
shows a detailed Q-dependence arising from changes in S(Q) (eq. |2J). The sharp maxima 
are Bragg-reflections of crystalline cubic ice as it is presented in detail in reference Jf^. 

The temperature dependence of the elastic intensity I(T,u = 0) is reported in the 
right panel of Fig. H Plotted data are taken at Q = 1.7 A -1 (Fig. EJa) and Q = 2.4 A -1 
(Fig. HJb). Figure EJc shows the signal integrated over the accessed Q range. It is evident 
from Figure 0] c that despite the intensity gain at the Bragg reflections (Fig. 0] a) and 
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intensity loss between (Fig. 0] b) the total scattering power in the accessed Q range does 
not show any anomalies in the temperature range of the supposed glass transition. In 
conclusion, the behaviour of I(T, to = 0) in Fig. 0] gives evidence of a redistribution of the 
signal within the elastic channel only, which is compatible with a recrystallisation of the 
LDA matrix. This recrystallisation is obviously slow at T < 150 K but very efficient above 



150 K 



T4L 0], causing the rapid gain in I{u = 0) shown in Fig. El The very onset of the 



recrystallisation matches with the supposed glass transition temperature T g « 135 K that 
is indicated by the dashed vertical line in Figure |U 



IV. b. Debye— Waller factor and mean— square displacement 



As it is outlined in section III the Debye-Waller factor W(Q,T) of a harmonic system 
obeys the simple relation W(Q,T) oc — T • Q 2 . Figures E] and El show the T- and the 
Q 2 -dependence of some selected W(Q,T) data sets for the samples SI and S2 measured 
on IN13. Since, the zero point oscillation of the scatterers is of no importance for the 
discussion of W(Q,T) here, the presented data sets are normalised to the signals measured 
at 2 K. It is obvious from Fig. E] that all phases, HDA, LDA, and ice I c , display a Q 2 - 
dependence and are (meeting the requirement of sufficiently high temperature) proportional 
to T. The linearity in T is equally fullfilled for the elastic intensity of S2 plotted in Fig. El 
up to temperatures close to the transformations of the phases. Grey shaded areas in Fig. El 
indicate the transformation regions, in which the structural correlations change, as it is 
shown in Fig. |U The Q 2 -dependence is as well fulfilled for HDA and LDA in S2 as it is for 
HDA, LDA and I c in SI, which is demonstrated by fits plotted as solid lines in Figs. Eland 
6 

According to eq. El the slopes of the fits correspond to the mean-square displacement 
{u 2 (T)) of the scatterers. They are shown for both SI and S2 in Fig. Despite the differ- 
ent levels of deuteration of SI and S2 the (u 2 (T)) of both samples in the HDA and LDA 
states are comparable. Please note that the onset of structural changes, which the samples 
encounter during the phase transitions, is demonstarted by the strong enhancement of sta- 
tistical uncertaintity. The comparable (u 2 (T)) values of both samples in the corresponding 
phases can be understood when taking into consideration that at the exploited temperatures 
(u 2 (T)) is strongly influenced by the mean velocity of sound (eq. EJ)- It has been shown that 
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FIG. 5: Debye- Waller factor W(Q, T) of HDA, LDA and ice I c measured with the partially pro- 
tonated sample SI on IN13. Left panel gives the temperature dependence of W(Q,T), right panel 
gives the Q-dependence of W(Q, T) plotted vs. Q 2 to stress the behaviour expected for a harmonic 
solid. The selected Q-values and approximate temperatures valid for all presented data are given 
in the top subfigures. Solid lines correspond to Q 2 -fits to the data whose slope is interpreted as 
the mean-square displacement (u 2 (T)) in Fig. [7J 

for the crystalline hexagonal ice Ih the level of deuteration does not change markedly the 
elastic constants or thermodynamic observables being determined by the vibrational proper- 



ties 



0,0, Ei. 



This means that the velocity of sound does not alter strongly with the level of 



deuteration. For example, a difference of only 3% was found in specific heat measurements 



of crystalline ice at T < 70 K as an effect of deuteration 
expected in the phases studied here. 



771 ] . A similar behaviour can be 
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FIG. 6: Debye- Waller factor W(Q, T) of HDA and LDA measured with the deuterated sample 
S2 on IN13. Left panel gives the temperature dependence of W(Q,T), right panel gives the Q- 
dependence of W(Q, T) plotted vs. Q 2 to stress the linear behaviour expected for a harmonic 
solid. The selected Q-values and approximate temperatures valid for all shown data are given in 
the top subfigures. Solid lines correspond to Q 2 -fits to the data whose slope is interpreted as the 
mean-square displacement (u 2 (T)) in Fig. [Tj Please note that original data have been grouped by 
having summed up two data points each. 

Assuming that the Debye model is applicable a Debye temperature Td can be extracted 
from fits to (u 2 (T)) via eq. ^2 Fit results are plotted as solid (SI) and dotted (S2) lines in 
Fig. [7| and the obtained Tb are given in Table |TI] The upper fitting limits are given by the 
transition temperatures which are obtained with sample S2 and shown in Fig. |UJ 

Another observable in neutron scattering from which (u 2 (T)) can be calculated is the 
density of states G{ui) (eg. EJ). Figure |H1 reports G{ui) measured on the time-of-flight spec- 



trometer IN6 



44, 



45l . 46]. The plotted G(tu) are calculated from data measured at 80 K 
(HDA), 127 K (LDA) and 160 K (ice I c ). Within the statistical accuracy of the measure- 
ments performed at different T no daviations from the presented G(u) could have been 
established. A self consistent multi-phonon calculation was applied to the data in addition 
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FIG. 7: Mean-square displacement (u 2 (T)) of HDA, LDA and ice I c determined from data taken 
with sample SI (■) and with sample S2 (□) shown in Figs. |S] and |BJ respectively. Solid lines 
correspond to fits with the Debye-model to data taken on SI, dotted lines correspond to fits to 
data on S2. Dashed lines are calculated with eq. |l]from the generalised density of states given in 
Fig. |HJ Please note that data of sample S2 have been shifted for a better visibility. 
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TABLE II: Debye-temperatures obtained from fits to the data shown in Fig. [7| and estimated 
from velocity of sound given in [&I . 



phase sample SI sample S2 from [83] 

HDA 230 ± 6 K 223 ±9K 230 K 

LDA 217±2K 202 ±3 K 203 K 

ice I c 224 ± 7 K 203 K 



to the standard correction procedures 



44. 



Since the data are obtained with a fully deuterated, coherently scattering sample the 
determined G{uj) is an approximation to the true density of states [3,0, IsjJ. To account 
for the perturbations by the coherent scattering a correction factor a = 0.7 is introduced. 
This factor is estimated independently comparing G{uS) of hexagonal ice Ih and ice XII 
measured on fully deuterated and fully protonated samples [82]. The Debye levels of the 
samples are indicatec 
experiments on ice I 



in Fig. |H1 and compared to the Debye levels obtained in ultrasonic 
83^ and Brillouin light scattering on ice Ih |84], both techniques giving 
comparable values. Taking a and G(u) from Fig. |H1 into account, a \ Uq^(T)J is calculated 
to the values presented as dashed lines in Fig. 

Within the statistics of the measurements and the accuracy of the applied data analysis 
the results obtained from backscattering and time-of-flight measurements fully agree. For 
example, in the framework of the Debye model both techniques indicate the highest average 
velocity of sound and, thus, the highest Td for HDA. Debye temperatures estimated via 
eq. HOI from results based on supersonic measurements |83] are included in Table HI and 
their respective Debye levels for G(u) are indicated in Fig. |H1 Although obtained with a 
different technique the Trj> match nicely the values observed in our elastic neutron scattering 
experiments giving also a higher value for Td for HDA 85^ . 

We would like to recall the fact that the coherent character of the sample introduces 
uncertainties into G(uS) although these are well controlled S^ . Thus, the excellent agreement 
of \Uq^(T)J with (u 2 (T)) calculated from the Debye-Waller factor W(Q,T) should not 
be over-interpreted. However, qualitatively the results on W(Q,T) and (u 2 (T)), which are 
determined by studying the elastic signal, and, on the other hand, the data on \Uq^(T)J, 
which is obtained from the inelastic response, give a coherent picture of the properties of 
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FIG. 8: Top: Generalised density of states G{oo) of HDA (■), LDA (□) and ice I c (♦) determined 
on the time-of-fnght spectrometer IN6 on a fully deuterated sample. All G{uj) are normalised to 
12 modes in the energy range fkv < 40 meV. Bottom: G(lu)/lu 2 emphasises the low energy range 
where a constant value is expected for harmonic solids. Please note that data at Huj < 1.5 meV are 
perturbed by the resolution of the spectrometer and therefore suppressed in the figure. Equally 
shown are the Debye-levels determined independently for fully deuterated (solid line) and fully 
protonated (dashed line) ice lh, from which a factor of 0.7 accounting for the coherency effect is 



estimated. Dotted line represents Debye-levels calculated from [83| for ice I n and from 
lh obtained by ultrasonic and light scattering experiments, respectively. 



841 ] for ice 
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HDA, LDA and ice I c . They mark them unequivocally as harmonic systems. 

Finally, we discuss the Debye temperatures determined by specific heat measurements 
[nl as 288 K for HDA, 305 K for LDA and 325 K for ice Ih- The fundamental difference is 
that Td is determined here for T <C Td, i.e. within the validity limits of the Debye model. 
For this range specific heat measurements on ice Ih give Td values of about 200 K with a 
difference of less than 3 % at T < 70 K between deuterated and protonated samples 7?j], 
as it was already indicated above. At elevated temparatures the specific heat reflects the 
details of the phonon density of states G(uS) beyond the acoustic region and the extracted 
Td becomes temperature dependent [72|. Taking only translational modes, i.e. phonons, 
into consideration the Td of ice Ih reaches a maximum at T pa 100 K with Td pa 300 K \T7\. 
Consequently, the basic difference of about 100 K between the T D presented here and cited 
in is due to the applied experimental techniques and the thermodynamic conditions of 
the measurements, but the data do agree in the general concept of harmonic solids 64 . t3| . 

The low value of Td of HDA given in can be, on the one hand, understood by 
the temperature sensitivity of Td in the specific heat measurements. At about 100 K not 



44, 



It is well established by 



only phonons but also librations contribute to C P (T) 
experiments that the librational band shifts toward lower energies the higher the density 
of the ice phase is. The shift of the librational band in HDA raises the specific heat 
and consequently pulls Td below the value of crystalline ice 4^, ^6. Ig^l . This effect is 
well comparable with the issue of deuterating a crystalline sample 



On the other 

hand, it must be stressed that the sample referred to as HDA in [lq was annealed before 
the measurements raising the question of the sample state. In the case of a mixture of 
HDA-type and LDA-type phases or even a contribution of ice XII j^] , 67,13] the density 
of states would decrease Td more efficiently when increasing T. 



VI. c. Excess of modes at low T 



So far we have not discussed the zero point oscillation of the scatterers in HDA, LDA 
and ice I c which is indicated in Fig. El Indeed, bearing in mind that the elastic intensity 
is coupled to the inelastic response (eqs. El El IH El) there is an inconsistency concerning 
the behaviour of HDA. Since, among all phases the lowest Debye-level is detected for HDA 
(Fig. El), it strictly requires the highest elastic intensity I(uj = 0) at T pa K. However, 
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HDA exhibits the lowest I {to = 0) measured in all our experiments (Fig. EJ). It must be 
notified that this inconsistency with I(u = 0) holds also 



"or ^results from other experiments 
Moreover, the Debye-Waller factors W(Q,T), the mean-square displacements (u 2 (T)), 



exploiting the dynamic response of the phases 



43]|4L 
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and the Debye temperatures Tb calculated from the elastic scan measurements itself display 
a dependence on temperature which is in full agreement with the G(u) determined by 
time-of-flight techniques. Consequently, the temperature evolution of the elastic intensity 
I(u = 0, T) of HDA contradicts its own value at very low temperatures I(u = 0, T — > 0). 

This self-contradicting behaviour of HDA can be resolved, if not only phonons, i.e. 
excitations following the Bose statistics, are inherent to the dynamic response but also 
non-Bose modes predominant at low temperatures are present. Such modes have been 
detected for the first time by three axis neutron scattering (TAS) (3, 54^ confirmed later 
and in terp reted as the contribution of two-level systems (TLS) to the inelastic response of 



HDA 



iterpretea as tn 

00,00 



. Remarkably, no TLS are observed in the LDA state. Please note, 



;hat the level of deuteration of the sample has no significant influence on the TLS intensity 

3,0,0- 



V. SUMMARY AND CONCLUSIONS 



We utilised the neutron back-scattering and time-of-flight techniques to gain informa- 
tion on the temperature dependence of the high-density (HDA) and low-density (LDA) 
amorphous states of ice and the crystalline phase I c . Particular focus was put on the phase 
transformations HDA — > LDA, LDA — > I c and the supposed glass-transition of LDA with 
T g w 135 K, the mean-square displacements of the scatterers in the stable states HDA, LDA 
and I c , and their relative zero point oscillations. 

Within the best energy resolution of 1 /ieV of the experiments no relaxation processes, 
giving evidence of a higher mobility of the scatterers, can be identified during the phase 
transformations. From results on coherently scattering samples, we may conclude that at or 
prior to T « 135 K, i.e., at the very T g discussed in literature, the recrystallisation of LDA 
into I c sets in - thus preventing any data analysis in terms of Mode Coupling Theory as it 
is oulined in jssj ] . 

The temperature and Q-dependence of the Debye-Waller factors and the resulting mean- 
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square displacements mark HDA, LDA and ice I c as harmonic solids. In the framework of 
the Debye model HDA reveals a higher Debye temperature, a higher average velocity of 
sound, and thus a lower phonon density of states at low energy when compared with LDA 
and ice I c . LDA and ice I c show rather similar dynamic properties. The conclusions drawn 
from the elastic scan data on the temperature dependence of all states are in full agreement 
with results obtained by other experimental techniques. 

The relative elastic intensities determined at low temperatures support the existence 
of additional modes in HDA not obeying the Bose statistics, i.e. this excess should not 
be confused with a Boson-peak as it is identified in glassy systems [41]. These modes have 



been observed in inelastic neutron scattering 
systems from light absorption experiments 



y, y, koi 



^confirmed and interpreted as two-level 



In the case of the protonated samples, the scattering from hydrogen is by a factor of 
~ 40 stronger than the scattering from oxygen. Consequently, the signal measured with 
samples S3 and S4 is not only determined by the dynamic properties of the water molecules 
but equally reflects the characteristics of the hydrogen sublattice. The absence of anomalies 
in the Debye-Waller factor excludes therefore any strong proton relaxation within the 
matrix. However, migration of lattice defects, as discussed in references 
detectable due to their too low concentration. 
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